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We apply a linearized stability analysis to simplified models of accretion-phase neutrino 
fluxes streaming from a supernova. We compare the results with recent numerical studies 
and find excellent agreement. This provides confidence that a linearized stability analysis 
can be further applied to more realistic models. 

1 Introduction 

Neutrino-neutrino interactions cause the neutrino flux evolution close to a supernova (SN) core 
to be nonlinear and numerically very challenging jTj. The flavor instability causing collective 
flavor conversions can be suppressed by the "multi-angle matter effect" [2]. This point was 
recently investigated numerically for an accretion-phase model where the matter density near 
the neutrino sphere is large, using a schematic description of the neutrino fluxes [5]. On the other 
hand, the flavor stability can also be investigated with a linearized stability analysis, avoiding 
an explicit solution of the equations of motion [3]. We apply this method to the models of 
Ref. [3] and find excellent agreement of the stable regime identified with either method. 

2 Linearized stability analysis 

We describe the neutrino flavor evolution in terms of matrices ^E.u,r where the diagonal ele- 
ments are the usual total number fluxes and the off-diagonal elements encode phase informa- 
tion [21 [5] • We label the angular dependence with u, in close relation with the neutrino emission 
angle -dn at the inner boundary radius R, u = sin^ -dfi = (1 — cos'^ i9r)r^ /R"^- For semi- isotropic 
emission at a "neutrino sphere" with radius R, the flux is uniformly distributed on < u < 1. 
The equations of motion are idr^E,u,r — [^E,u,r,^E,u,r], with the Hamiltonian 

H^,... = f ^ + ^/2 GpN,) ^ + ^ d^.' dE' ( -J— l) 0^, , 

\2-C; / Vu,r 47rH Jq J_^ \Vu,rVu',r J 

where is the neutrino mass-squared matrix, the matrix of net chargcd-lepton densities 
which in the flavor basis is Hf = diag(ne— ng, Ufi—Ufi, nr—rif) and Vu.r is the radial projection of 
neutrino velocity at the radius r. Antineutrinos are represented through negative-energy modes 
(E < 0) and negative negative fluxes in the matrices 'i>E,u,r- This sign convention simplifies 
the formalism and obviates any distinction between neutrinos and antineutrinos. 
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Henceforth we drop the exphcit subscript r to denote the r-dependence of all quantities. In 
the two flavor scenario one can write: 

Tr (i>E.u , ^E.u.B ^ ^E.u.R c 
'^E,u — ^ 1 ^ o_B,ti , 

where -F^'^ are the differential neutrino fluxes at the inner boundary radius R for the e and x 
flavors. The flux summed over all flavors, Tr 'i>E,m drops out of the equations of motion and is 
conserved in our free-streaming limit. The "swapping matrix" 

Se.u Se.u 

encodes the flavor evolution with initial conditions s = 1 and = 0. 

We expand the Hamiltonian for large distances from the core and small mixing angle. After 
dropping its trace we flnd 

Hvac _ M!„-i^±^ ("co^Se sin20\ 1 ^/l OW 
nE,u 2E 2 \sm2e -cos2e) " 2 \0 -l) \ 2 7-2 

He,u V2Gf Q Qji'„ ^-(^Q -ij 1,^+2 7^ 

where A = ^/2Gf {ne—rie). We write the neutrino-neutrino part concisely in the form H'^'^^ = 
fir Jq <iu' ^{u + u') dE' gE^uSE'.u', where fir = fJ-nR'^/^r^ encodes the strength of the 
neutrino-neutrino interaction with the parameter fijj = y^GpiF^" — i^]^'")/47ri?'*. We further 
define the dimensionless flavor difference spectrum gE.u = {F% „ ^ — -Ff „ ^) / {F'^' — F'^ ) with 
the nornicLlizcition in tlic antincutrino sector J_^dE ^'^Qe.u ~ 

— 1. The integration over 

neutrinos (positive energies) gives dE du gE,u = (Fj^' — F'^'')/{F'^ — F'^) = 1 + e, with 
e being asymmetry of the spectra. 

Next we expand the equations in the small-amplitude limit jS*! <C 1 which implies, to linear 
order, s = 1. After switching to the variable uj = /S.m^ /2E for the energy modes one finds [4] 

'ydrSuj,u = + u{\ + e^)] Suj.u - fi du' dcu' {u + u') g^j'u' S'w'.u' • 



Here A = XR^/2r^ encodes the imprint of multi-angle matter effect. Except for the additional 
two powers of this quantity describes the SN density profile and scales approximately as 
fir oc r~^. 

Writing solutions of the linear differential equation in the form Suj,u = Qi^.u e^'^'' with 
complex frequency — ^ + \k and eigenvector Qu},u leads to the eigenvalue equation [4], 

(u) + uX - fl) Quj,u = fi j du' du' (u + u') g^iu' Qlj'm' , 

where A = \ + £fi. The solution has to be of the form Qi^,u = (A + Bu) /{uj + uX — Vl). Solutions 
exist if fi^^ = /i ± ^/IqI2, where /„ = J duj du g^^u u" /{uj + uX — fl). The system is stable if all 
^l are purely real. A possible imaginary part, k, is the exponential growth rate. 



2 



HANSE 2011 



3 Results 



We aim at comparing the linearized stability analysis with the numerical solutions of Ref. [3] 
who numerically solved the neutrino flavor evolution for a 10.8 Mq model at various post bounce 
times. They confirmed the multi-angle matter suppression of self-induced flavor conversion, but 
also found partial conversions at a large radius for the models 200 ms < tpb ^ 300 ms. 

We use the same schematic half-isotropic and monochromatic spectra, leading to the simple 
form Qi^^u = —S{uj + ujq) + {I + s) 5(w — loq). The integrals /„ can now be evaluated analytically. 
Then it is easy to find a solution (7, k) for each pair (/i, A). Figure [T] shows the region where 
K 7^ for two snapshots together with the k isocontours. Wc also show the "SN trajectory" in 
the (/X, A) plane, i.e. essentially the density profile as a function of radius because oc r^^. 

Our results agree with the numerical solutions of Ref. [3] for all models. Whenever the 
numerical solutions find no flavor conversion, our SN trajectory indeed stays clear of the unstable 
regime. Conversely, when it briefly enters the unstable regime as in the left panel of Fig. [Tl we 
reproduce the onset radius for partial flavor conversion of Ref. [3] . The linear stability analysis 
correctly explains the numerical results. 
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Figure 1: Contours of k and the trajectory of SN (thick red line) at t — 300 ms (left) and 
400 ms (right) post bounce for a 10.8 A/q model discussed in Ref. [3]. 



It is interesting that in principle the SN 
trajectory can enter the instability region 
twice. As a toy model we consider the density 
profile A ~ 0.43 ^ with half-isotropic emission 
at R = 10 km and fir ^ 7 x 10^ km"^ R'^/2r'^. 
In Fig. [5] we show K{r) and the evolution of 
the off-diagonal element \S\. Indeed \S\ oscil- 
lates and grows in the unstable regime, only 
oscillates when k = 0, and then grows again 
during the second instability crossing. It re- 
mains to be seen if there are realistic density 
profiles where such a multiple instability sit- 
uation exists in practice. 
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Figure 2: Growth rate k and off-diagonal ele- 
ment \S\ for a toy model (see text). 
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4 Conclusions 



The nonlinear neutrino flavor evolution in the SN environment can be a challenging numerical 
task even when it only consists of post-processing the output of a self-consistent SN simulation, 
not to mention solving self-consistently the multi- flavor neutrino transport. The latter is not 
necessary if collective oscillations do not happen in the critical region below the shock wave. 
The question if a given SN model with concomitant neutrino fluxes is stable against self-induced 
flavor conversion can be answered with a linearized stability analysis [4]. Of course, if the model 
is unstable, one needs to solve the equations numerically to find the final outcome. However, 
since neutrino fiuxes during the accretion phase may well be stable because of the multi-angle 
matter effect, a linearized fiavor stability analysis is here a useful tool. 

We have applied this method to the models studied in Ref. [3] and compared with the 
outcome of their numerical solutions. The results are very encouraging in that we can perfectly 
account for the results of the numerical simulations and can also predict the onset radius for 
those cases where partial flavor conversion occurs at a large radius. 

Meanwhile we have applied this method to an accretion-phase model with realistic energy 
and angle distributions [5]. We find that a stability diagram in the form of our Fig. [l]is an 
excellent tool to summarize the flavor stability situation of SN neutrino fiuxes. 
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